Crossover from Selberg's type to Ruelle's type zeta function in classical kinetics. 
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The decay rates of the density-density correlation function are computed for a chaotic billiard 
with some amount of disorder inside. In the case of the clean system the rates are zeros of Ruelle's 
Zeta function and in the limit of strong disorder they are roots of Selberg's Zeta function. We 
constructed the interpolation formula between two limiting Zeta functions by analogy with the case 
of the integrable billiards. The almost clean limit is discussed in some detail. 
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It is natural to assunie j-that chaotic billiards with a 
small amount of disordeiui3 are good models for ballistic 
cavities which have been employed in a number of recent 
experiments, see Ref. ^. Such a model is interesting be- 
cause the disorder (in two dimensions) can be character- 
ized by one parameter: the elastic scattering time r. The 
mixing properties of this model in two important limits 
T ^ and r — > OO were known in the literature, see Refs. 
o,@. In the present work we discuss the crossover from 
one limit to the other for some two-dimensional billiard. 
The case of the three-dimensional billiards is more com- 
plicated. Particularly, the uniform scattering in three di- 
mensions leads to very fast resonant mixing. At the end 
of the paper we provide the generalization of our results 
for the case of three dimensions. 

Let us focus attention on the eigenvalues and the cigen- 
modes of the kinetic equation for the distribution func- 
tion /(r, (p) of non- interacting particles inside a two- 
dimensional billiard. This function is defined on the 
constant energy manifold \v^\ = v =const, and v^ — 
{v cos((/()), V sin(0)). The precise form of the kinetic equa- 
tion depends on the details of the impurity potential, but 
we are going to investigate two models 
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where /(r) = / /(r, (p) d0/(27r). The above equation has 
to be solved with mirror boundary conditions f{r, cj)) = 
f{r, 2a{r) — tt — 0), where r is taken on the boundary 

of the billiard and n = I cos(a), sin(Q;) ) is normal to the 

boundary. Equation {M has a special solution fo{f, (f>) — 
const for all values of r and we will ignore it in the rest 
of the paper. 

In both models, the collision integral conserves energy. 
The first model corresponds to uniform scattering in all 
directions and the second model is valid if small angle 
scattering is dominant. Let us look for solutions propor- 
tional to e"''"*. The eigenvalues s„ of the kinetic equation 



are so-called mixing rates, or decay rates of the density- 
density correlation function when t — > or Ruelle's res- 
onances when r ^ OO. 

These resonances can be found as zeros of the spectral 
determinant Z(s). Let us start to compute|^(s) in the 
limit of pure chaos. Cvitanovic and EckardtEl have com- 
puted the spectral determinant for the axiom A system, 
but the result is the expansion over the unstable periodic 
orbits and it seems to be valid for wide class of systems. 
Therefore in the limit t — > oo the spectral determinant 
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where v is the velocity. This expression contains the sum 
over the primitive periodic orbits p taken with repetition 
r. In the case of billiards, the action for the periodic 
trajectory is just the length Ij^r. Each oscillating term in 
the sum in Eq. (||) is weighted by the stability amplitude 
which behaves on average as 
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where the first equality defines Apr, and A is the Lya- 
punov exponent of the billiard. The Zeta function given 
by Eq. (|) is of Ruelle's type. 

When the disorder is strong, the kinetic equation can 
be transformed into the diffusion equation for / 
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which has to be solved with boundary conditions nV/ = 
0. Thjs equation allows one to find the decay of modes 
with / ^ 0. The decay rates for modes with / = for all 
r should be computed in a different way. 

In the limit r ^ we can use the "semiclassical" ap- 
proximation for equation 



(fc2 + v2)/ = 
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where fc^ = 2s/{v'^t). The logarithm of the Selberg's 
type Zeta function is again the sum over periodic orbitaj'L 
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Here the phase factors correspond to the case when 
Eq. (1q) should be solved with the Neumann boundary 
condition and Maslov's indexes vanishJlj The natural 
question to ask is whether it is possible to compute the 
decay rates of the modes with / ^ for all values of r 
by constructing a suitable Zeta function? 

We can understand the connection between different 
types of Zeta functions by making use of the following 
approximation 
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where we have neglected the repetitions of periodic orbits 
and fluctuations of the Lyapunov exponent in Eqs. (j^) 
and (||). Let us denote the eigenvalues of the wavenum- 
ber in Eq. (||) as g„, n = 1 . . . oo. We therefore suggest 
that eigenvalues of the kinetic equation are 
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FIG. 1. The decay rates of the density-density correlation 
function move on the complex plane when disorder decreases, 
as shown by arrows. The termination point is Res ~ A/2, 
where A is Lyapunov exponent. 

Let's examine first the integrable case. Model 1 f©r 
the square billiard was_solved by Atland and GefentI, 
and Agam and Fishmann, who modeled the short range 
random potential by random spheres or circles. For 
the square billiard of size L the spatial dependence of 
the density /(r) is ^exp(igr), where q is such that 
sm(qxL) sm{qyL) — 0, and the sum is over four possi- 
ble directions of q. Then the values of q are "quantized" , 
and we will denote them q„, and the modes with / ^ 
can be numbered 
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where A is the mean Lyapunov exponent. In order to ob- 
tain the second expression we have noted that the func- 
tion A{x) should have zeros a;„ = ig„ — ^ and we have 
made use of Eq. (g). We would like to emphasize that 
Eq. (nG) has nothing to do with the convergence prop- 
erties of the Zeta functions. It based on the fact that 
Ruelle's type Zeta function Eq. (g) and Selberg's type 
Zeta function are almost the same function A{x) taken 
with different arguments. The second case in Eq. (|l(]| ) 
means that there is a shift of A/2 between the zeros of 
Ruelle's type Zeta function and the zeros of Selberg's 
type Zeta function. Such a shift was observed in the 
problem of quantum and classical scattering in a three 
disk problem, compare Figs. 2.14 and 3.6 of Ref. O and 
see Refs. |2|,|3| for details. The main idea of this work is 
to obtain known Zeta functions as the two limits of one 
kinetic Zeta function. 
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where the sum is over four possible directions of qn 
tegration over leads to the equation for s„ 
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and the corresponding Zeta function is 



log(Z(s))=> J-^ 
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and 



where p is not a single orbit but the resonant toriEil 
the connection between s and q is as in Eq. (E2[). 

In the case of model 2, the solutions are still propor- 
tional to e**'', but the angular dependence is different. 
The solution with / ^ is the "ground state" of 
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/n = 0, 
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because the real parts of the decay rates are positive. 
Surprisingly, Eq. (|l2| ) gives a numerically good approx- 
imation for s„ for this model. Other "angular" modes, 
which have / = 0, are very different for models 1 and 2. 

It is not easy to compute eigenmodcs of Eq. (nl) for 
the integrable billiards which have other than rectan- 
gular shapes. For such cases Eq. ( |l3| ) becomes an in- 
terpolation formula for the kinetic Zeta function. One 
should only replace the pre-expanjential factor by the am- 
plitude from the Berry - TaborEj formula. For example, 
the resonant tori for the circular billiard of radius R are 
numbered by the winding number M and by the num- 
ber of vertices n, have length Imti — nLMn/T^, where 
Lmu — 2TTRsm{7rM/n). Then one can use Eq. (|lj) after 
the replacement p — > Mn, and L -^ Li\jn, see Rc f. [iq . 

Combining together Eqs. (||), (||), (p^, and (|13|) we 
can introduce the kinetic Zeta function as 
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p r—1 



where 



' pr 



r, Xprr/2 < 1 

Z/Apr-i AprT / -Z ^ J- 



(15a) 



(15b) 



For At < 1/2 the kinetic Zeta function coincides with 
Selberg's type Zeta function Eq. (g) in the domain of the 
complex s plane \st\ <C 1. For Ar > 1/2 the kinetic Zeta 
function becomes independent of r and coincides with 
Ruelle's type Zeta function Eq. (0) in the domain of the 
complex s plane \st\ ^ 1. 

The interpolation formula Eqs. (nS) for the kinetic Zeta 
function implies the following interpolation formula for 
the decay rates 
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where g„ are the eigenvalues of the wavenumber in 
Eq. (H). There is a gap between the last two expressions 
Sn|r=2/A-o - s„|r=2/A+o ^ ^^ / ilnv) , which is numeri- 
cally small for most cases. The motion of the decay rates 
on the complex plane is schematically shown in Fig. P. 
In the limit of strong disorder some of the s„ are on the 
real axis, and the imaginary part of s„ becomes non-zero 
when TQnV = 1. Then s„ moves along the arc and stops 
when T = A/2. 

Equations (|o|) and (|6|) show that a chaotic system is 
qualitatively different from a diffusive system from the 
point of view of the position of Ruelle's resonances s„ on 
the complex plane. In the chaotic limit all resonances 
lie on a line parallel to the imaginary axis. The disorder 



induces motion of the resonances toward the real axis as 
was found by Agam and Fishman.cl 

The interpolation formula between Ruelle's type and 
Selberg's type Zeta exists only if the diffusion modes 
transform to the so-called Frobenius - Perron modes as 
the disorder goes to zero. This has not yet been proven 
for our case. The main difficulty is that the diffusion 
modes are selected from all kinetic modes by the condi- 
tion / ^ 0. At the same time Frobenius-Perron modes 
are selected by the choice of the functional space. How- 
ever, in other systems one can consider the diffusion 
modes as modes of the Frobenius - Perron operator .113 

Some additional information might be obtained from 
the properties of the propagator of Eq. (|l|), which can be 
written as a path integral for model 2 



G(f, (/), fo,(/)o,i) = / L'[V'](5(f-fo - / v^dt) 
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where ip{Q) = </>o, t}}{tj ~ Q) + ijj{tj + 0) = 2aj, ... 
, ip{t) — (/), and the path tq + J^ v^dt touches the 
boundaries n times at the points ri, . . . , Vn, at the times 
ti, . . . , t„. The angle aj is the direction of the tangent to 
the boundary at the reflection point Tj . The trace of the 
propagator Eq. (nw known also as the return probability 
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where ip{t) = ip{0) = (f> and J^ ip'^dt is defined as in 
Eq. (0). 

The propagator Eq. (|l^) should interpolate between 
the Frobenius - Perron operator in the limit t —f 00 and 
the diffusion operator in the limit r — > 0. Then the trace 
of this propagator Eq. (fsl) should provide us with a sys- 
tematic way to compute the interpolation formula for the 
Zeta function, because — log Z{s) = /^ e"* t~^p{t)dt. 
Here the sign of st in the Laplace transform is positive, 
because we want the roots of the Zeta function to have 
the meaning of the decay rates. 

In the limit of weak disorder r — > 00, one may hope 
to obtain the small corrections oc ^ to the Frobenius - 
Perron operator, and therefore to Eq. (g). Particularly 
one may expect to obtain the additional "stabilization" 
of the periodic orbits through the disorder. Let's con- 
sider the vicinity of the periodic orbit p in phase space. 
The path in such a vicinity can be described by the coor- 
dinate x{t) — vt along the orbit, by the coordinate y(t) 
normal to the orbit, and by the deviation of the direction 
of motion 4){t). The position of the particle at the end of 
the path and at the beginning of the path are connected 
by 
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where the orbit p crosses the billiard Up times. In other 
words the orbit consists of Up segments of length Lpj. 
When the particle is going along the segment j, it can 
be scattered by the disorder at small angle 9j, and then 
the rest of the path is distorted too. The cumulative 
change of the end of the path is given by the sum in 



the right hand side of Eq. (|19D, where Mpj is the mon- 
odromy matrix of the piece of the orbit consisting of the 
segments L 



pj + l: ■ 
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One can see immediately from 
Eq. ( |19| ) , that the stability amplitude of the closed path 
y{t) = 2/(0), (j){t) = 0(0) is independent of 6*1, . . . , ^n^, 
and therefore it is independent of r. Therefore, there 
are no l/r corrections to the Zeta function Eq. (0) and 
our interpolation formula Eq. (|I5|) is independent of r for 
T > 2/A. 

In the case of the three dimensional billiards the ef- 
fect of the disorder is different because the scattering 
becomes three-dimensional and the distribution function 
depends on the three coordinates and two angles. For 
the cubic billiard the spatial dependence of the density 
is again ^e*'''', where the sum is over the six orthogo- 
nal orientations of q, and the modes are selected by the 
rule siii{qxL) sin(gj,i) shi{qzL) — 0. Then, the dispersion 
relation for the analog of model 1, (uniform scattering), 
can be found in Ref. [iq : 

tan(qwT) 

where qvr < tt. In other words there are no modes with 
qvT > TT, and / ^ 0, where the bar means the average 
over the solid angle. Equation ( |2C| ) describes the diffu- 
sion modes for small r, but it cannot be used for large r. 
If the mode has q close to tt/{vt) then the decay of such 
a mode is very fast s oc (tt — qvT)~^. 

The model with small angle scattering in three dimen- 
sions is the Fokker - Planck equation for the distribution 
function, which should be solved together with mirror 
boundary conditions on the billiard boundary. The so- 
lutions inside the cubic billiard have the same dispersion 
s{q) as in the case of the square billiard, if / ^ 0. There- 
fore one may hope that Eq. (|l^) gives the interpolation 
of the Zeta function of the Fokker - Planck equation for 
modes with / ^ 0. 

In summary, we have constructed the interpolation for- 
mula for the Zeta function of the kinetic equation, in both 
"chaotic", Eqs. ([15|), and "integrable" , Eq. (|l^) cases. 
From the mathematical point of view our kinetic Zeta 
function interpolates between Ruelle's and Selberg's Zeta 
functions. Our formulas are independent of the partic- 
ular choice of the collision integral for two-dimensional 
billiards and are suitable for small angle scattering in 
three dimensions. 
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